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Abstract 

We analyze the exact behavior of the renormahzation group flow in one-dimensional 
clock-models which undergo first order phase transitions by the presence of complex in- 
teractions. The flow, defined by decimation, is shown to be single-valued and continuous 
throughout its domain of definition, which contains the transition points. This fact is in 
disagreement with a recently proposed scenario for first order phase transitions claiming 
the existence of discontinuities of the renormahzation group. The results are in partial 
agreement with the standard scenario. However in the vicinity of some fixed points of the 
critical surface the renormalized measure does not correspond to a renormalized Hamil- 
tonian for some choices of renormahzation blocks. These pathologies although similar to 
Griffiths-Pearce pathologies have a different physical origin: the complex character of the 
interactions. We elucidate the dynamical reason for such a pathological behavior: entire 
regions of coupling constants blow up under the renormahzation group transformation. 
The fiows provide non-perturbative patterns for the renormahzation group behavior of 
electric conductivities in the quantum Hall effect. 
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The behavior of the renormahzation group flow around first order transition points 
has been a controversial matter for years. The conventional scenario is based on a smooth 
behavior of the renormahzation group flow and the existence of a discontinuity fixed point 
whose attraction domain contains the transition surface and has relevant exponents of the 
form y = D [1] [2] , D being the dimensionality of the system. The singularities associated 
with first order transitions are generated by infinite iterations of renormahzation group 
transformations in the thermodynamic limit. 

However, recently, some authors claimed to have numerical evidence of the failure of 
the conventional picture [3] [4]. They proposed a new scenario where the renormahzation 
group fiow is discontinuous at the first order transition surface. The ambiguity arises be- 
cause the renormahzation group transformation associates to any system on the transition 
surface as many different renormalized Gibbs measures as phases coexist at the transition. 

More recently, van Enter, Fernandez and Sokal [5] rigorously proved that the second 
picture turns to be false. Their result states that for systems with bounded fluctuating 
variables and absolutely summable real Hamiltonians, the renormahzation group transfor- 
mation, when properly deflned by some blocking procedure, is single-valued and continuous 
in a domain of the space of parameters which includes flrst order transitions points. The 
only pathology which arises for some systems (e.g. Ising model in D > 2) is that for some 
renormahzation group prescriptions (majority rule, decimation, etc.) the renormahzation 
group transformation is not defined at all in a neighborhood of first order transition points, 
as previously pointed out by Griffiths-Pearce [6] and Israel [7]. The numerical results of 
ref. [3] might be understood as artifacts of the truncation of the local renormalized Hamil- 
tonians used in the Monte Carlo renormahzation group analysis [8]. The method used in 
ref. [4] does not rely on truncation of local interactions but assumes the existence of only 
one relevant perturbation at the fixed point [9] and in case of Ising model there are two: 
one associated to the temperature T = 1/ {3 and another to an external magnetic field H 
[2] [8]. 

In this note we analyze those issues in some exactly soluble models recently introduced 
by two of us [10]. The simplicity of the models allow us to analyze exact renormahzation 
group flows in flnite-dimensional coupling spaces. Although they exhibit very peculiar 
properties due to the presence of complex interactions, we can extract some rigorous lessons 
about the behavior of the renormahzation group . The physical interest of those models 
comes from the relevance of their (cr-model) higher dimensional generalizations in statistical 
mechanics and quantum fleld theory (cr-models with similar complex couplings arise e.g. 
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in the eflFective description of the quantum Hall effect [11]). They can also be considered as 
one-dimensional versions of Chiral Potts models which have been a recent focus of interest 
[12]. 

The complex character of the interactions of the models makes possible the existence 
of first order phase transitions though they have only short range interactions (in fact 
nearest- neighbor), which is impossible for one-dimensional systems with real Hamiltonians 
(see ref. [13] for a recent review). Another unusual property of the models is that the 
correlation length at the (first-order) transition surface is infinite which implies a very 
peculiar behavior of finite size effects in the vicinity of the critical points [14] . 

Our results for the three and four-state clock models are in agreement with the con- 
ventional scenario for the renormalization group near first order critical points. The renor- 
malization group transformation is single-valued and continuous throughout its domain 
of definition which contains transition points, and has discontinuity fixed points with the 
peculiarity that the correlation length at those fixed points is infinite (^ = oo). 

However, there are some pathologies of the systems at the points of the coupling space 
where the renormalization group transformation is not defined. They arise because of the 
appearance of renormalized measures without renormalized Hamiltonians and by no means 
are related to the ambiguities of the discontinuous scenario. In these sense they are analo- 
gous to Griffiths- Pearce singularities. In the present case there is also a dynamical reason 
behind such a pathological behavior: the blow up of the corresponding renormalization 
group trajectories at infinite under a finite number of (fractional) iterations. The renor- 
malization group flow generates a (local) one-parameter group of local transformations [15] 
but there is no global renormalization group. If we exclude the pathological region of pa- 
rameters in the vicinity of first order transition points the renormalization group is globally 
defined. This standard dynamical behavior is the reason for such a pathological behavior. 
However, the physical origin of these pathologies is different from that of Griffiths-Pearce 
singularities. In this case is the existence of complex interactions interactions which makes 
possible the appearance of renormalized measures which vanish for a very large sets of 
configurations and cannot be described in terms of renormalized Hamiltonians. In the 
present models the pathological values of the couplings are non-generic and non universal. 
They depend on the renormalization group prescription and might be thought as artifacts 
of the choice of the prescription. However, we find pathologies for any prescription based 
on decimation. 
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The renormalization of systems with complex topological ^-terms (cr-models, gauge 
theories, etc) has been always problematic from the numerical point of view, and semi- 
classical considerations based in instanton calculus cannot be taken for grant. The exact 
results obtained in our models shed some light in the possible (pathological) behavior of 
the renormalization group for those systems and provide a non-perturbative pattern which 
seems to pick up the essential behavior of ^-terms under renormalization. 

We consider the g-state clock of a classical spin variable s„ fluctuating among the 
g-roots of unity, 

Sn = (cos— — ,sm— — ) pn = 0,...,q-l (1) 
in one-dimensional space with interacting Hamiltonian [10] 

N 

f3H=-^ {JSn ■ Sn+1 - ie- {Sn X S^+l)) 
n=l 

= - f JCOS {pn - Pn+l) + iesin {pn - Pn+l) 

n=l V ^ ^ 



(2) 



where /3 denotes as usual the inverse of the temperature. The eigenvalues of the transfer 
matrix T 

q~l 

Xk = ^ exp {Jcos(27rn/q') + z£sin(27rn/g) — ik{27rn/q)} (3) 

n=0 

cross each other an infinity number of times. The crossings of leading levels correspond 
to first order transition points with infinite correlation length [10]. In the case of 3-state 
clock model the three eigenvalues of the transfer matrix 

Ao = e-^ + 2e~^ cos(6'), Ai = e-^ -2e"^ cos(6' + 7r/3), A2 = e^^ -2e"^ cos(6'-7r/3) (4) 

have leading level crossings at the transition temperatures 9 = (2m -|- l)7r/3, where 
9 = £\/3/2. Single alternate decimation induces the following renormalization group trans- 
formation 

^, 2, e'^'^ + 2e--^ ^, 2ei sm9 - e'-^ sm29 

J = - log — -, 9 = arctan — j . (5) 

3 (e~2'^ -h 4e'^ -h 4e~ cos36')2 2e2 cos^ -I- e~'^cos26' 

There is an additional free energy renormalization which is not relevant for the analysis of 
the critical behavior of correlation functions. The transformation (5) has to be understood 
as the projection of the renormalization group transformation on the plane of coupling 
constants (J, 9). 
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The flow is univocaUy defined everywhere and exhibits a very interesting dynamical 
behavior. There are not traces of the ambiguities associated to discontinuous scenarios. 

The critical surfaces 9 = (2m + l)7r/3 are invariant under renormalization group transfor- 
mations as well as the subcritical ones 9 — 2m7r/3. There are three kinds of fixed points 
(see Fig. 1) 



oo 


J = 


0,9 = 2n7r/3 


(stable) 


infinite temperature 


-^0 


J = 


oo 


(unstable) 


zero temperature 


pn 
^ c 


J = 


|log2,6' = (2n + l)7r/3 


(saddle point) 


transition temperature 



The third class of fixed points are sitting on the transition curves and attract all 
transition points. There is one relevant direction at those points which is tangent to the 
renormalized trajectory 

q3J/2 ^ 

9 = arccos(— — ) + (2m + 1)^ (6) 

and flows towards the inflnite temperature fixed points . The correlation length vanishes 
at the infinite temperature fixed point F^ and becomes infinite at the points F" and Fq 
of the critical surface. 

The only second order critical points where the continuum limit can be obtained by 
scaling are the zero temperature fixed points Fq. The existence of a line of fixed points 
leads to the existence of different quantum mechanical systems [16]. 

Another interesting aspect of the renormalization group arc the critical exponents. 
The linearized renormalization group tranformations around the different fixed points are 

^^^ = (o o) ^^" = (o 2) ^n^ = (o i)- 

Therefore, the infinite temperature fixed points F^ are stable and have two irrelevant 
directions (A^ = 2^^ = 0, i = 1, 2). 

The second order transition points at zero temperature have two marginal perturba- 
tions (y = 0), one along the line of fixed points J = and another tangent to a renormalized 
trajectory fiowing away towards the other types of fixed points. 

The fixed points on the first order transition surface F" are saddle points with one 
relevant direction (Ai = 2^^ = 2) flowing along the renormalized trajectory (6) towards the 
the infinite temperature fixed points F^ and F^^ and one irrelevant direction (A2 = 2^^ = 
0) along the critical line 9 = {2n + l)7r/3. The existence of one relevant perturbation with 



y = 1 agrees with the picture advocated by the standard scenario for discontinuity fixed 
points of first order transitions except for the fact that the correlation length is infinite. In 
general, the number of relevant perturbations with critical dimension equal to the space 
dimension D = 1 has be equal to the number of phases minus one. In this case it implies 
that the critical exponent u = 1. 

The only amazing behavior of the renormalization group flow arises for 

^<^o(6') = ^log2cos(6'-(2m + l)|). (8) 

where some eigenvalues of the transfer matrix T become negative. One single iteration of 
the renormalization group transformation maps any point of this region into a point of the 
region J > Jo{0) where the transfer matrix is positive (T > 0), and the renormalization 
group transformation is singular at the transition points J = — |log2,^ = (2n + 1)tv/3 
which are mapped by (5) into J = oo,6 = (2n + l)7r/3. Since our decimation procedure 
transforms the transfer matrix T into T^, which is always positive the corresponding 
renormalization group transformation maps the points of the region J < Jo{0) into those 
of the region J > Jo{0). The renormalized trajectory (6) is the borderline between those 
domains. 

The singularity at the points J — oo,6 = {2n + 1)tt/3 is due to the diagonal character 
of T^. The vanishing of non-diagonal entries implies that the renormalized measure is not 
Gibbsian and does not correspond to any regular renormalized Hamiltonian. 

In order to analyze a more natural flow for Hamiltonians with T ^ we deflne a 
renormalization group transformation by double decimation (see Fig. 2), which leaves 
invariant the domain where T is not non-positive. The corresponding renormalization 
group transformation 

2 , e^-^ + 6 + 2e-3-^/2 cos 39 

J = o log 



3 3(4cosh^3J/2 + l + 4cosh3J/2cos36')^ 

,, 2 cosh 3 J/2 sin ^ — sin 2^ 

6 = arctan — 

2 cosh 3 J/2 cos 6* + cos 26* 

exhibits a similar global {universal) behavior for the T > region J > Jo{d) (see Fig. 
1), and it is not singular at J = — |log2,^ = (2n + l)7r/3. In fact, these points become 
flxed (unstable) points. However, the transformation (9) is again singular at the points 
satisfying that 

e^-^ + 6 + 2e-3-^/^ cos 3^ = 0, (10) 
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which are mapped into the infinite fine J = — oo. 

The analysis of the renormahzation group flow for e^"^ + 6 + 2e~^"^/^ cos 3^ < requires 
an analytic extension of the model to the complex plane of the coupling constant J [17]. 
The only really pathological points are those of the curve (10) which can be characterized 
by the fact that the diagonal elements of vanish. In such a case the renormalized 
measure vanishes for a large set of clock variables configurations and cannot be described 
in terms of Gibbsian positive weights associated to a renormalized Hamiltonian. 

Therefore, the pathological behavior of systems without non-negative transfer matrix 
can be improved but not completely cured by this change of the blocking procedure. 
However, we remark that the set of pathological Hamiltonians is non-generic in the space 
of couplings constants: it has codimension one or two depending on the blocking procedure. 

The universal behavior for ,/ > Jo{6) also holds for the continuous fiow associated 
to the transformation T ^ which is only valid for T > 0. The integration of the 
corresponding differential equations 

f 2 , , 2,^ 2e-3>^/\ • e'^/2 + 2e-'^cos6',, , , ^ . ^ x 

J= -e ■^AologAo + -(l + — )cot^^ ^ ^ . ^ (AilogAi - AalogAa) 

6 6 COS U 3^3 sm t) 

a/3 

d= 4^e^/2(3cos^(AilogAi - AalogAa) - \/3sin6'(2Ao log Aq - Ai log Ai - A2 log A2)) 

lo 

(11) 

yields a continuous 

fiow with the same global characteristics and the same picture for first order transition 
points (see Fig. 1). 

The behavior of the renormahzation group is essentially changed by the presence of 
a external real magnetic field H. When H ^ there are not first order transition points 
and Lee- Yang theorem holds in this case. The renormahzation group can be exactly 
analyzed by introducing a new coupling which breaks the internal rotation symmetry of 
the clock model. The only first order critical points Fc belong to the surface H = and a 
new relevant direction appears at along the renormalized trajectory associated to the 
external magnetic perturbation. 

In the case of the 4-state model whose critical behavior was analyzed in [10], the 
renormahzation group transformation induced by decimation is not well defined unless 
we introduce a new interaction term of the form Ji cos7r(pn, — p^+i) in the Hamiltonian. 
Therefore we need to consider a new family of models with Hamiltonians 

N 

Pn' = - XI ("^ f -Pn+l) + JlCOSn{pn-Pn+l)+i£Sm^{pn-Pn+l)^ ■ (12) 
n=l 
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The eigenvalues of the transfer matrix are given by 

Ao = 2e'^i cosh J + 2e~''^ cos e Ai = 2e'^' sinh J + 2e~'^^ sin e 

(13) 

A2 = 2e'^i cosh J - 2e~"^i cos e A3 = 2e'^i sinh J - 2e~"^i sin e 

and crossings of leading eigenvalues correspond again to first order transition points with 
infinite correlation length. 

An alternate single decimation defines the following renormalization group transfor- 
mation 

, 1, l + e*"^icosh2J 

J = - log -7-, 

2 ^ c^'^i + cos 2e 
,_1 (l + e-4-^icos2£)(l + e4'^icosh2J) r^^^ 
^ ~4 2cosh2J + 2cos2e 

, 1 1 + cosh2Jcos2£ 

£ = - arccos — — — . 

2 cosh 2 J + cos 2e 

In this case a new problem arises. The transformation (14) is not defined for couplings of 

the domain 

B2 = {{J, Ji, s); e^-^i < - cos 2s, e > (Black Hole) (15) 

in terms of renormalized hamiltonians of the type (12) with real couplings J, Ji, e. However, 
there exist an analytic extension of Ti' to the complex planes of J and Ji that preserves 
the hermitian character of the transfer matrix; and the renormalized hamiltonians of the 
points of B2 belongs to such a class of hamiltonians [17]. Thus, the only real pathology 
arises for couplings in the border of the black hole 

dB2 = {(J, Ji,e);e^-^i = -cos2£, e > ^}. 

Those pathological points are non-generic and vary with the definition of the renormal- 
izacion group transformation. A common characteristic of those points is that they yield 
renormalized measures which vanish for a large set of clock variables configurations and 
cannot be described in terms of a Gibbs measure associated to a renormalized Hamiltonian. 

Apart from such a pathological zone the renormalization group is well defined ev- 
erywhere and none ambiguity arises at first order transition surfaces. Because of the 
periodicity in s it suffices to consider the region < £ < ^. The critical surfaces in that 
region are 

g-J+2Ji ^ _V2cos(£ + t) t < £ < ^ (16a) 

^4 4 2 

J<2Ji, £=| (166) 



and the renormalization group fixed points are listed below (see also Fig. 3) 
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(saddle point) 
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F* 


J = 
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(saddle point) 


zero temperature 
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Ft 
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Once again positivity provides a useful information about the renormalization group 
flow. The domain where the transfer matrix is positive (T > 0) is located below the surface 



Jo(Ji, £) = max{arcosh (e ^ cos £) , arcsinh (e ^sin£)}. (17) 

The points above such a surface J < Jo(Ji,£) where T is not non-positive are mapped 
by one single iteration of the renormalization group transformation into points below the 
surface, J > Jo(Ji,£), by the same reason that in the case of the 3-state clock model. A 
double iteration would cure this anomaly but then the domain of definition of the renormal- 
ization group transformation in terms of Hamiltonians of the type (12) with real couplings 
J, Ji,e is smaller than that of a single alternate decimation (14). The renormalization 
group flow for the region where T is not non-positive J < Jq will be discussed in [17]. Let 
us concentrate here on the region with positive transfer matrix J > Jq. On the first order 
transition surface the only fixed points are F^, F* and Fj. Although the attraction do- 
main of all those points contains transition points, only F^ has one eigenvalue 2 {y = 1) of 
the linearized renormalization group transformation as required in the standard Nienhuis- 
Nauenberg picture for discontinuity fixed points. In our case it implies that the critical 
exponent v = 1, 

The renormalization group is well defined on the domain of the transition surface 
with J < Jq. Transition points with J < 2Ji in the plane £ = 7r/2 are attracted by FJ, 
transition points with J = 2 Ji move towards the fixed point F* along the renormalized 
trajectory defined by the edge J = 2Ji,£ = 7r/2 of the transition surfaces (16a) and 
(166), and the remainder points of the transition surface with J > 2Ji flow towards the 
flxed point F* . A similar behavior is observed for points below the transition surface with 
J > 2Ji, which reach the black hole domain B2 after a flnite number of iterations. This is 
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one of the reasons why the (complex) region (black hole) for T^"' grows with n. The 
new pathological points for T^"^ are points that already reached the region B2 in less than 
n iterations of . 

On the other hand, the continuous flow defined for systems with T > by the equa- 
tions 

j = -^(Ai log Ai + A3 log A3)e'^' cosh J - ^(Aq log Aq + A2 log A2)e'^i sinh J 

J ^ (AplogAp + A2logA2) (AplogAp - A2logA2) ^ tang ^ tanhJ 
^~ Se-^icoshJ 8e--^icos£ ^2 2 ^ ' 

e = — (cos£(Ai log Ai - AslogAa) - sin£(AplogAp - A2logA2)) 

is defined everywhere below the surface (17). In particular, it is well defined in the patho- 
logical points of dB2. The only pathological property of the points of the domain B2 is 
that they have been swallowed up by the line J = 00, Ji = — oo,£ e (7r/4, 7r/2] in a finite 
time t < 2. This feature explains why there is a restriction on the domain of definition 
of T^. In fact, all the points below the transition surface will reach such a line in a finite 
time, which means that any given point of such a domain does not belong to the domain 
of for n large enough. The points at the transition surface (16a) with J > 2Ji, J > Jq 
are also attracted by the fixed point F* but in an infinite time (The intersection curve of 
the critical surface and J = Jp is the renormalized trajectory joining F* and (see Fig. 
3)). Therefore, the renormalization group transformation is well defined in such points for 
any kind of decimation blocking. 

The pathologies observed so far in the definition of the renormalization group trans- 
formation can be then understood as a pure consequence of the dynamical behavior of 
a renormalization group flow with a very strong attractive domain. For any continuous 
flow with such a strong domain attractor the corresponding (local) one-parameter group of 
local transformations cannot be implemented to a global group of transformations [15]. In 
the present case the transformation T* cannot be globally defined for any value of t 7^ 0. 
However, if we exclude the region below the transition surface (16a) the renormalization 
group fiow defines a global one-parameter group of transformations for any value of t e R. 

In summary, besides the above dynamical system explanation, we have shown three 
sources of pathological behavior of the renormalization group, all of them motivated by 
the complex nature of the interacting terms of the Hamiltonian: 
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i) The presence of negative eigenvalues in the transfer matrix yields an odd behaviour 
under decimation renormalization group transformations which generates pathologies. In 
particular, a continuum renormalization group cannot be implemented in that case. 

ii) In the 4-states model there are some values of the parameters (region B2) where 
the transfer matrix is positive but the renormalization group is not defined. This can 
happen because the correspondence between eigenvalues of and the parameters of the 
renormalized Hamiltonian is not analytic. 

iii) It is possible to extend the domain of definition of the renormalization group 
transformation by analytic continuation in the space of coupling parameters [17]. In this 
way the set of pathological Hamiltonians can be reduced to non-generic hypersurfaces with 
non-trivial codimension in the space of coupling parameters (e.g. the boundary of B2). 

In any case, although the Hamiltonians of these models are complex the Boltzmann 
weights can be paired to give rise to real contributions. The main difference with standard 
systems with real couplings is that some of the weights might vanish. However, the config- 
urations with vanishing contributions are non-generic and should not be confused with the 
pathological renormalized measures, where the configurations with vanishing contributions 
become generic (e.g. the transfer matrix has null all non-diagonal entries). 

One might conjecture that the pathological behavior of renormalization group trans- 
formation observed in higher dimensional systems is generated by similar dynamical prop- 
erties. In fact, the same behavior occurs in the renormalization group flow of the qf-state 
clock models with Hamiltonians 



which are the natural generalization of the systems (2). 

Finally, we remark that the renormalization group flow of the g = 3 model (Fig. 1) is 
similar to the expected flow of dissipative and Hall axy conductivities in the quantum 
Hall effect when described in terms of a cr-model [18]. Although the clock model can be 
considered as a one-dimensional discretization of such a cr-model, it was not expected to 
belong to the same renormalization class and provide an exact non-perturbative pattern 
for the renormalization group behavior of those conductivities (We thank Alexei Morozov 
for pointing out such a connection). 

The most remarkable property of this pattern is that in the thermondynamic limit 
the dissipative conductivity would vanish and the Hall conductivity would be quantized 




(19) 
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as it is observed in the quantum Hall effect (see fig.l of the first and third papers of Ref. 
[18]). These facts raise the conjecture that a solution of the strong CP-problem of QCD 
might be found in a similar way in terms of a renormalization flow of the 6 parameter 
towards quantized (unobservable) values 6 = 2nn. In the continuum limit all values of 
9 are consistent but in a given (perhaps fundamental) discretisation the topological CP- 
breaking couplings flow in the infrared to the quantized CP-preserving values 27rn. 
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Figure Captions 

Fig. 1.: Renormalization group flow for the 3-state model. The renormalization group 
transformatfon is universal for J > Jo(Ji, e) (T > 0). The renormalization group 
trajectories are integral curves of the continuous flow deflned by equation (11). 
Besides the zero temperature fixed points Fg, which are not displayed in the 
picture, there are two types of fixed points: infinite temperature fixed points 
F^{J = 0,9 = 2n7r/3), and discontinuity fixed points {J = |log2, ^ = 
{2n + l)7r/3). The later behave as attractors for the transition points on the lines 
9 = (2n + l)7r/3. In the domain J < Jo(Ji,e) where T < 0,the renormalization 
group is non-universal. In fact, it becomes singular at the points with J ~ 
-2/3 log 2, 9 ={2n + l)7r/3 and the curve e^'^ + 6 + 2e-^^/'^ cos 3^ = (dashed 
curve) for single and double decimations, respectively. 

Fig. 2.: Single and double decimations. 

Fig. 3.: Renormalization group flow for 4-state model. Shadowed surfaces correspond 
to transition points satisfying the positivity condition T > 0. For simplic- 
ity, we only show the spectrum of fixed points and renormalized trajectories 
of the renormalization group transformation. The boundary of the domain S2, 

dB2 = {( J, Ji, e), e^'^i = — cos 2£, £ > ^} where the renormalization group trans- 
formation is not defined can be seen on the left back corner of the (e, e-4"^i , 0) 
plane. 
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